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A Comparison of the One-Dimensional Bridge Hydraulic Routines from: 
HRC-RAS. HEC-2. and WSPRO 


INTRODUCTION 

The hydraulics of flow through bridges is an important aspect of computing water surface 
profiles. The computation of accurate water surface profiles through bridges is necessary in 
flood damage reduction studies, channel design and analysis, and stream stability and scour 
evaluations. There are several one-dimensional water surface profile computer programs 
available for performing this type of computation. The most widely used of these programs are: 
HEC-2 (HEC, 1991) and WSPRO (FHWA, 1990). The Hydrologic Engineering Center (HEC) 
has recently released a new program for computing one-dimensional water surface profiles, 
called HEC-RAS (HEC, 1995). The purpose of this study was to evaluate the effectiveness of 
the new bridge hydraulics routines in HEC-RAS at sites with extensive observed data, and to 
compare HEC-RAS to HEC-2 and WSPRO, with respect to bridge modeling performance. 

Detailed data for 22 bridge sites were obtained from the US Geological Survey (USGS, 1978). 
The USGS collected extensive data for bridge crossings over wide, densely vegetated flood 
plains between 1969 and 1974. The data has been published in the form of Hydrologic 
Investigation Atlases, which are available to the public. Of the 22 sites, 13 were used for this 
study. Several of the sites had to be omitted because of sparse water surface measurements in the 
vicinity of the bridge. A few other sites were omitted due to inadequate bridge geometry and 
layout information. Of the 13 remaining data sets, 17 events were analyzed; that is, 8 sites had 1 
flood each; and 4 sites had 2 floods each. Almost all of the events were class A low flow (open 
channel, subcritical flow through the bridge opening), while three of the events had water 
surfaces higher than the bridge low cord on the upstream side of the bridge. 


OVERVIEW OF BRIDGE HYDRAULICS 

In general, all three of the models compute water surface profiles by solving the energy equation, 
utilizing the standard step procedure. However, one of the main differences among these models 
is how they compute the water surface profile through a bridge. A description of the bridge 
hydraulics for each of the three models follows: 


HEC-RAS 

The HEC-RAS bridge routines utilize four user-defined cross sections in the computation of 
energy losses due to the structure. A plan view of the basic cross section layout is shown in 
Figure 1. 
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Cross Section Locations 


Cross section 1 is located sufficiently downstream from the structure so that the flow is not 
affected by the structure (i.e., the flow has fully expanded). This distance should generally be 
determined by field investigation during high flows. If field investigation is not possible, then 
there are two sets of criteria for locating the downstream section. The USGS has established a 
criterion for locating cross section 1 a distance downstream from the bridge equal to one times 
the bridge opening width (the distance between points B and C on Figure 1). Traditionally, the 
Corps of Engineers criterion has been to locate the downstream cross section about four times the 
average length of the side constriction caused by the structure abutments (the average of the 
distance from A to B and C to D in Figure 1). The expansion distance will vary depending upon 
the degree of constriction, the shape of the constriction, the magnitude of the flow, and the 
velocity of the flow. Both criteria should be used as rough guidance for placing cross section 1. 
The user should not allow the distance between cross section 1 and 2 to become so great that 
friction losses will not be adequately modeled. If the modeler thinks that the expansion reach 
will require a long distance, then intermediate cross sections should be placed within the 
expansion reach in order to adequately model friction losses. The ineffective flow option can be 
used to limit the effective flow area of the intermediate cross sections in the expansion reach. 


Approach Section 





y Exit Section 


Figure 1. HEC-RAS Cross Section Locations at a Bridge 
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Cross section 2 is located immediately downstream from the bridge (i.e., within a few feet). 

This cross section should represent the effective flow area just outside the bridge. 

Cross section 3 should be located just upstream from the bridge. The distance between cross 
section 3 and the bridge should be relatively short. This distance should only reflect the length 
required for the abrupt acceleration and contraction of the flow that occurs in the immediate area 
of the opening. Cross section 3 represents the effective flow area just upstream of the bridge. 
Both cross sections 2 and 3 will have ineffective flow areas to either side of the bridge opening 
during low-flow and pressure-flow profiles. In order to model only the effective flow areas at 
these two sections, the modeler should use the ineffective flow area option at both of these cross 
sections. 

Cross section 4 is an upstream cross section where the flow lines are approximately parallel and 
the cross section is fully effective. Because flow contractions can occur over a shorter distance 
than flow expansions, the distance between cross section 3 and 4 should be roughly one times the 
average width of the opening. However, this criterion for locating the upstream cross section 
may result in too short a reach length for situations where the width of the bridge opening is very 
small in comparison to the floodplain. An alternative criterion would be to locate the cross 
section a distance upstream equal to the average contraction width (the average of the distance 
from A to B and C to D in Figure 1). 

During the hydraulic computations, the program automatically formulates two additional cross 
sections inside of the bridge structure. The geometry inside of the bridge is a combination of the 
bounding cross sections (sections 2 and 3) and the bridge geometry. The bridge geometry 
consists of the bridge deck and roadway, sloping abutments if necessary, and any piers that may 
exist. The user can specify different bridge geometry for the upstream and downstream sides of 
the structure, if necessary. Cross section 2 and the structure information on the downstream side 
of the bridge are used as the geometry just inside the structure at the downstream end. Cross 
section 3 and the upstream structure information are used as the bridge geometry just inside the 
structure at the upstream end. 


Hydraulic Computations Through the Bridge 

The bridge routines in HEC-RAS allow the modeler to analyze a bridge with several different 
methods without changing the bridge geometry. The bridge routines have the ability to model 
low flow (Class A, B, and C), low flow and weir flow (with adjustments for submergence on the 
weir), pressure flow (orifice and sluice gate equations), pressure and weir flow, and highly 
submerged flows (the program will automatically switch to the energy equation when the flow 
over the road is highly submerged). 
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Low Flow Computations 


Low flow exists when the flow going through the bridge opening is open channel flow (water 
surface below the highest point on the low chord of the bridge opening). For low flow 
computations, the program first uses the momentum equation to identify the class of flow. This 
is accomplished by first calculating the momentum at critical depth inside the bridge at the 
upstream and downstream bridge sections. The cross section with the higher momentum 
(therefore most constricted section) will be the controlling section in the bridge. If the two 
sections are identical, the program selects the upstream bridge section as the controlling section. 
The momentum at critical depth in the controlling section is then compared to the momentum of 
the flow downstream of the bridge when performing a subcritical profile (upstream of the bridge 
for a supercritical profile). If the momentum downstream is greater than the critical depth 
momentum inside the bridge, the class of flow is considered to be completely subcritical (i.e., 
class A low flow). If the momentum downstream is less than the momentum at critical depth, in 
the controlling bridge section, then it is assumed that the constriction will cause the flow to pass 
through critical depth and a hydraulic jump will occur at some distance downstream (i.e., class B 
low flow). If the profile is completely supercritical through the bridge, then this is considered 
class C low flow. 

Class A low flow. Class A low flow exists when the water surface through the bridge is 
completely subcritical (i.e., above critical depth). Energy losses through the expansion (sections 
2 to 1) are calculated as friction losses and expansion losses. Friction losses are based on a 
weighted friction slope times a weighted reach length between sections 1 and 2. The average 
friction slope is based on one of the four available alternatives in the HEC-RAS, with the 
average-conveyance method being the default. This option is user selectable. The average length 
used in the calculation is based on a discharge-weighted reach length. Energy losses through the 
contraction (sections 3 to 4) are calculated as friction losses and contraction losses. Friction and 
contraction losses between sections 3 and 4 are calculated in the same way as friction and 
expansion losses between sections 1 and 2. 

There are three methods available for computing losses through the bridge (sections 2 to 3): 

- Energy Equation (standard step method) 

- Momentum Balance 

- Yamell Equation 

The user can select any or all of these methods in the computations. If more than one method is 
selected, the user must choose either a single method as the final solution or direct the program to 
use the method that computes the greatest energy loss through the bridge as the answer at section 
3. This allows the modeler to compare the answers from several techniques all in a single 
execution of the program. Minimal results are available for all the methods computed, but 
detailed results are available for the method that is selected as the final answer. A detailed 
discussion of each method follows. 
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Energy Equation (standard step method) : 


The energy based method treats a bridge in the same manner as a natural river cross section, 
except the area of the bridge below the water surface is subtracted from the total area, and the 
wetted perimeter is increased where the water is in contact with the bridge structure. As 
described previously, the program formulates two cross sections inside the bridge by combining 
the ground information of sections 2 and 3 with the bridge geometry. As shown in Figure 2, for 
the purposes of discussion, these cross sections will be referred to as sections BD (Bridge 
Downstream) and BU (Bridge Upstream). 


(,3) (BU) (Bd (2) 




Figure 2. Cross Sections Near and Inside the Bridge 


The sequence of calculations starts with a standard step calculation from just downstream of the 
bridge (section 2) to just inside of the bridge (section BD) at the downstream end. The program 
then performs a standard step through the bridge (from section BD to section BU). The last 
calculation is to step out of the bridge (from section BU to section 3). 
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Momentum Balance Method: 


The momentum method is based on performing a momentum balance from cross section 2 to 
cross section 3. The momentum balance is performed in three steps. The first step is to perform 
a momentum balance from cross section 2 to cross section BD inside the bridge. The equation 
for this momentum balance is as follows: 


y4 V + 
BD 


^bdQbd 


gA 


BD 


^ 2^2 - ^ 2 y ,2 " 


Ml 

gA. 




( 1 ) 


where: A^, 

= Active flow area at section 2 and BD, respectively 

^p2 

= Obstructed area of the pier on the downstream side 

Y2, Ybd 

= Vertical distance from water surface to center of gravity of 

flow area A 2 and Ago, respectively 


= Vertical distance from water surface to center of gravity of 

wetted pier area on downstream side 

P2» PbD 

= Velocity weighting coefficients for momentum 

Ql’ Qbd 

= Discharge 

g 

= Gravitational acceleration 

Ff 

= Force due to friction 

K 

= Force due to weight of water in the direction of flow 


The second step is a momentum balance from section BD to BU (see Figure 2). The equation for 
this step is as follows: 
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The final step is a momentum balance from section BU to section 3 (see Figure 2). The equation 
for this step is as follows: 
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where: = Drag coefficient for flow going around the piers. 

The momentum method provides detailed output for the cross sections inside the bridge (BU and 
BD) as well as outside the bridge (2 and 3). The user has the option of turning the friction and 
weight force components off; the default is to include these two terms in the momentum balance. 


Yamell Equation : 

The Yamell equation is an empirical equation that is used to predict the change in water surface 
from just downstream of the bridge (section 2 of Figure 2) to just upstream of the bridge (section 
3). The equation is based on approximately 2600 lab experiments in which the researchers 
varied the shape of the piers, the width, the length, the angle, and the flow rate. The Yamell 
equation is as follows (Yamell, 1934): 

H,,^2K(K+ 10(0 - 0.6)(a + 15a^) — (4) 


where: F/ 3.2 
K 
(0 

a 

F 2 


The drop in water surface from section 3 to 2. 

Yamell's pier shape coefficient 

Ratio of velocity head to depth at section 2 

Obstmcted area of the piers divided by the total unobstmcted area. 

Velocity downstream at section 2 


The computed upstream water surface elevation is simply the downstream water surface 
elevation plus Hj.j. With the upstream water surface known, the program computes the 
corresponding velocity head and energy elevation for the upstream section (section 3). When the 
Yamell method is used, hydraulic information is only provided at cross sections 2 and 3 (no 
information is provided for sections BU and BD). 

The Yamell equation is sensitive to the pier shape (K coefficient), the pier obstmcted area, and 
the velocity of the water. The method is not sensitive to the shape of the bridge opening, the 



shape of the abutments, or the width of the bridge. Because of these limitations, the Yarnell 
method should only be used at bridges where the majority of the energy losses are associated 
with the piers. 

Class B low flow. Class B low flow can exist for either subcritical or supercritical profiles. For 
either profile, class B flow occurs when the profile passes through critical depth in the bridge 
constriction. For a subcritical profile, the momentum equation or the energy equation can be 
used to compute an upstream water surface (section 3 of Figure 2) above critical depth and a 
downstream water surface (section 2) below critical depth. For a supercritical profile, the 
bridge is acting as a control and is causing the upstream water surface elevation to be above 
critical depth. Momentum or energy can be used to calculate an upstream water surface above 
critical depth and a downstream water surface below critical depth. If the user is running a mixed 
flow regime profile, the program will proceed with backwater calculations upstream, and later 
with forewater calculations downstream from the bridge. 

Whenever class B flow is found to exist, the user should run the program in a mixed flow regime 
mode. The mixed flow regime mode is capable of calculating a subcritical profile upstream from 
the bridge and a supercritical profile downstream from the bridge. Also, any hydraulic jumps can 
be located, if they exist. 

Class C low flow. Class C low flow exists when the water surface through the bridge is 
completely supercritical. The program can use either the energy equation or the momentum 
equation to compute the water surface through the bridge for this class of flow. 


High Flow Computations 

The HEC-RAS program has the ability to compute high flows (flows that come into contact with 
the maximum low chord of the bridge deck) by either the Energy equation (standard step method) 
or by using separate hydraulic equations for pressure and/or weir flow. The two methodologies 
are explained below. 

Energy Equation /standard step method! : 

The Energy based method is applied to high flows in the same manner as it is applied to low 
flows. Computations are based on balancing the energy equation in three steps through the 
bridge. Energy losses are based on friction and contraction and expansion losses. Output from 
this method is available at the cross sections inside the bridge as well as outside. 

As mentioned previously, friction losses are based on the use on Manning's equation. 

Contraction and expansion losses are based on a coefficient multiplied by the absolute change in 
velocity head. The Energy based method performs all computations as though they are open 
channel flow. At the cross sections inside the bridge, the area obstructed by the bridge piers. 
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abutments, and deck is subtracted from the flow area and additional wetted perimeter is added. 
Occasionally the resulting water surfaces inside the bridge (at sections BU and BD) can be 
computed at elevations that would be inside of the bridge deck. The water surfaces inside of the 
bridge reflect the hydraulic grade line elevations. 


Pressure and Weir Flow Method : 

A second approach for the computation of high flows is to utilize separate hydraulic equations to 
compute the flow as pressure and/or weir flow. The two types of flow are presented below. 

Pressure Flow Computations. Pressure flow occurs when the flow comes into contact with the 
low chord of the bridge. Once the flow comes into contact with the upstream side of the bridge, 
a backwater occurs and orifice flow is established. The program will handle two cases of orifice 
flow; the first is when only the upstream side of the bridge is in contact with the water; and the 
second is when the bridge constriction is flowing completely full. The HEC-RAS program will 
automatically select the appropriate equation, depending upon the flow situation. For the first 
case (see Figure 3), a sluice gate type of equation is used (FHWA, 1978): 


Q = C, [2g (T3 


Z 

2 




2 ^ 


)]^ 


(5) 


where: Q 

Q 

■^BU 

Z 


Total discharge through the bridge opening 
Coefficient of discharge for pressure flow 
Net area of the bridge opening at section BU 
Hydraulic depth at section 3 

Vertical distance from maximum bridge low chord to the mean 
river bed elevation at section BU 


The discharge coefficient C^, can vary depending upon the depth of water upstream. Values for 
Cj range from 0.35 to 0.5, with a typical value of 0.5 commonly used in practice. The user can 
enter a fixed value for this coefficient or the program will compute one based on the amount that 
the inlet is submerged. A diagram relating to Y 3 /Z is shown in Figure 4. 
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Figure 3 Example of a bridge under sluice gate type of pressure flow 


As shown in Figure 4, the limiting value of Y 3 /Z is 1.1. There is a transition zone somewhere 
between Y3/Z = 1.0 and 1.1 where free surface flow changes to orifice flow. The type of flow in 
this range is unpredictable, and equation 5 is not applicable. 

In the second case, when both the upstream and downstream side of the bridge are submerged, 
the standard full flowing orifice equation is used (see Figure 5). This equation is as follows: 


Q = CAs^ 


( 6 ) 


where: C = Coefficient of discharge for fully submerged pressure flow. 

Typical value of C is 0.8. 

H = The difference between the energy gradient elevation upstream and 

the water surface elevation downstream. 


A - Net area of the bridge opening. 


I 

I 
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Figure 5. Example of a bridge under fully submerged pressure 




Typical values for the discharge coefficient C range from 0.7 to 0.9, with a value of 0.8 
commonly used for most bridges. The user must enter a value for C whenever the pressure flow 
method is selected. The discharge coefficient C can be related to the total loss coefficient, which 
comes from the form of the orifice equation that is used in the HEC-2 computer program (HEC,, 
1991): 



2gH 




K 


(7) 


where: K = Total loss coefficient 
The conversion from to C is as follows: 


C = 



( 8 ) 


The program will begin checking for the possibility of pressure flow when the computed tow 
flow energy grade line is above the maximum low chord elevation at the upstream side of the 
bridge. The user has the option to tell the program to use the water surface, instead of energy, to 
trigger the pressure flow calculation. Once pressure flow is computed, the pressure flow answer 
is compared to the low flow answer, the higher of the two is used. 


Weir Flow Computations. Flow over the bridge, and the roadway approaching the bridge, is 
calculated using the standard weir equation (see Figure 6): 


Q = CLH^'^ 


(9) 


where: Q = Total flow over the weir 

C = Coefficient of discharge for weir flow 


L = Effective length of the weir 

H = Difference between energy upstream and road crest 
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Figure 6. Example bridge with pressure and weir flow 

The approach velocity is included by using the energy grade line elevation in lieu of the upstream 
water surface elevation for computing the head, H. 

For high tailwater elevations, the program will automatically reduce the amount of weir flow to 
account for submergence on the weir. This is accomplished by reducing the weir coefficient 
based on the amount of submergence. Submergence corrections are based on a trapezoidal weir 
shape or optionally an ogee spillway shape. The total weir flow is computed by subdividing the 
weir crest into segments, computing L, H, a submergence correction, and a Q for each section, 
then summing the incremental discharges. The submergence correction for a trapezoidal weir 
shape is from "Hydraulics of Bridge Waterways" (Bradley, 1978). Figure 7 shows the 
relationship between the percentage of submergence and the flow reduction factor. The percent 
of submergence is computed as the depth of the tailwater above the weir divided by the upstream 
energy head above the weir, times 100. 
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Discharge Reduction for Submerged Flow 



Figure 7. Factor for reducing weir flow for submergence 


When the weir becomes highly submerged the program will automatically switch to calculating 
the upstream water surface by the energy equation (standard step backwater) instead of using the 
pressure and weir flow equations. The criterion for switching to energy based calculations is user 
controllable. A default maximum submergence is set to 0.95 (95 percent). 

Combination Flow 

Sometimes combinations of low flow or pressure flow occur with weir flow. In these cases, an 
iterative procedure is used to determine the amount of each type of flow. The program continues 
to iterate until both the low flow method (or pressure flow) and the weir flow method have the 
same energy (within a specified tolerance) upstream of the bridge (section 3). The combination 
of low flow and weir flow can be computed with any of the low-flow methods mentioned 
previously, except the momentum based method. 
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HEC-2 


HEC-2 computes energy losses caused by structures such as bridges in two parts. One part 
consists of the losses that occur in reaches immediately upstream and downstream from the 
bridge where contraction and expansion of the flow is taking place. The second part consists of 
losses at the structure itself and is calculated with either the normal bridge method or the special 
bridge method. Cross section locations are the same as described previously for HEC-RAS. 

Losses due to contraction and expansion of flow between cross sections are determined by 
standard-step profile calculations. Manning's equation is used to calculate friction losses, and all 
other losses are described in terms of a coefficient multiplied by the absolute value of the change 
in velocity head between adjacent cross sections. 


Normal Bridge Method 

The normal bridge method handles a bridge cross section in the same manner as a natural river 
cross section, except that the area of the bridge below the water surface is subtracted from the 
total area, and the wetted perimeter is increased where the water is in contact with the bridge 
structure. The user is required to enter the two cross sections inside the bridge. The bridge deck 
is described either by entering the constant elevations of the top of roadway and low chord, or by 
specifying a table of roadway stations and elevations, and corresponding low chord elevations. 
Pier losses are accounted for by the loss of area and the increased wetted perimeter of the piers, 
as described in terms of cross section coordinates. 


Special Bridge Method 

The special bridge method computes losses through the structure for either low flow, pressure 
flow, weir flow, or for a combination of these. The profile through the bridge is calculated using 
hydraulic formulas to determine the change in energy and water surface elevation through the 
bridge. 

Low Flow. The procedure used for low flow calculations in the special bridge method depends 
on whether the bridge has piers. Without piers, the low flow solution is accomplished by 
standard-step calculations as in the normal bridge method. The transfer to the normal bridge 
method is necessary because the equations used in the special bridge method for low flow are 
based on the obstruction width due to the piers. Without piers, the special bridge solution would 
indicate that no losses would occur. 

For a bridge with piers, the program goes through a momentum balance for cross sections just 
outside and inside the bridge to determine the class of flow. The momentum calculations are 
handled by employing the following momentum relations based on the equations proposed by 
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Koch-Carstanjen [Eichert/Peters, 1970] [Koch-Carstanjen, 1962]. 


/«, - m^, + 




A, - -^A 


pi 


= 


8^2 


m3 - mp3 + 




( 10 ) 


where: A„ A 3 = flow areas at upstream and downstream sections, 

respectively 

Aj = flow area (gross area - area of piers) at a section within 

constricted reach 


Ap„ Ap3 


Yn y2. Ys 


m„ m 2 , m 3 

mp„ mp3 

Cd 


Ypn Yp3 


Q 

g 


obstructed areas at upstream and downstream sections, 
respectively 

vertical distance from water surface to center of gravity of 
A„ Aj, A 3 , respectively 

A,yi, h -^2 A 3 y 3 , respectively 

Apjp, and Ap 3 yp 3 , respectively 

drag coefficient equal to 2.0 for square pier ends and 1.33 
for piers with semicircular ends 

vertical distance from water surface to center of gravity of 
Ap, and Ap 3 , respectively 

discharge 

gravitational acceleration 


The three parts of the momentum equation represent the total momentum flux in the constriction 
expressed in terms of the charmel properties and flow depths upstream, within and downstream 
of the constricted section. If each part of this equation is plotted as a function of the water depth, 
three curves are obtained representing the total momentum flux in the constriction for various 
depths at each location. The desired solutions (water depths) are then readily available for any 
class of flow. The momentum equation is based on a trapezoidal section and therefore requires a 
trapezoidal approximation of the bridge opening. 

Class A low flow occurs when the water surface through the bridge is above critical depth, i.e., 
subcritical flow. The special bridge method uses the Yamell equation for this class of flow to 
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determine the change in water surface elevation through the bridge. As in the momentum 
calculations, a trapezoidal approximation of the bridge opening is used to determine the areas. 


V 

= 2KiK + lOw - 0.6) (a + ISa'*) — (11) 

2g 


where: H 3 = 

drop in water surface from upstream to downstream sides of the 
bridge 

K 

pier shape coefficient 

0 ) = 

ratio of velocity head to depth downstream from the bridge 
obstructed area 

a = 

total unobstructed area 

V3 = 

velocity downstream from the bridge 


The computed upstream water surface elevation is simply the downstream water surface 
elevation plus H 3 . With the upstream water surface elevation known, the program computes the 
corresponding velocity head and energy elevation for the upstream section. 

Class B low flow can exist for either a subcritical or supercritical profile. For either profile, 
class B low flow occurs when the profile passes through critical depth in the bridge constriction. 
For a subcritical profile, critical depth is determined in the bridge, a new downstream depth 
(below critical) and the upstream depth (above critical) are calculated by finding the depths 
whose corresponding momentum fluxes equal the momentum flux in the bridge for critical 
depth. The program does not provide the location of the hydraulic jump. A supercritical profile 
eould be computed starting at the downstream section with a water surface elevation X. For a 
supercritical profile, the bridge is acting as a control and is causing the upstream water surface 
elevation to be above critical depth. Momentum equations are again used to recompute an 
upstream water surface elevation (above critical) and a downstream elevation below critical 
depth. 

Class C low flow is computed for a supercritical profile where the water surface profile stays 
supercritical through the bridge constriction. The downstream depth and the depth in the bridge 
are computed by the momentum equations based on the momentum flux in the constriction and 
the upstream depth. 

Pressure Flow. The pressure flow computations use the orifice flow equation of U.S. Army 
Engineer Manual 1110-2-1602, "Hydraulic Design of Reservoir Outlet Structures," [USAGE, 
1963]. 
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where: H 

— 

difference between the energy gradient elevation upstream and 
tailwater elevation downstream 

K 


total loss coefficient 

A 

= 

net area of the orifice 

g 

= 

gravitational acceleration 

Q 

= 

total orifice flow 


The total loss coefficient K, for determining losses between the cross sections immediately 
upstream and downstream from the bridge, is equal to 1.0 plus the sum of loss coefficients for 
intake, intermediate piers, friction, and other minor losses. The section on loss coefficients 
provides values for the total loss coefficient and shows the derivation of the equation and the 
definition of the loss coefficient. 

Weir Flow. Flow over the bridge and the roadway approaching the bridge is calculated using the 
standard weir equation: 

Q = CLH^^ (13) 


where: C 


coefficient of discharge 

L 

= 

effective length of weir controlling flow 

H 

= 

difference between the energy grade line elevation and the roadway 
crest elevation 

Q 

= 

total flow over the weir 


The approach velocity is included by using the energy grade line elevation in lieu of the upstream 
water surface elevation for computing the head, H. Where submergence by tailwater exists, the 
coefficient 'C’ is reduced by the program [Bradley, 1978]. Submergence corrections are based on 
a trapezoidal weir shape or optionally an ogee spillway shape. A total weir flow, Q, is computed 
by subdividing the weir crest into segments, computing L, H, a submergence correction and Q for 
each segment, and summing the incremental discharges. 
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Combination Flow. Sometimes combinations of low flow or pressure flow occur with weir 
flow. In these cases a trial and error procedure is used, with the equations just described, to 
determine the amount of each type of flow. The procedure consists of assuming energy 
elevations and computing the total discharge until the computed discharge equals, within 1 
percent, the discharge desired. 


WSPRO 

WSPRO computes the water surface profile through a bridge by solving the energy equation. 
Cross sections are located as shown in Figure 8. 


Flow 


2: Bridge 
Opening 


B 


1: Approach 
Section 






I 3F: Full 
jLi Valley 


I 3; Bridge 
I Opening 


4: Exit Section 

/ 


B 


Figure 8. Cross section locations for WSPRO bridge hydraulics 


Computations of the water surface profile requires the user to enter a minimum of four cross 
sections. These include sections 4, 3F, 3, and 1. Sections 4, 3F, and 1 are unconstricted full 
valley sections, while section 3 is the bridge opening section. Cross section 2 represents an 
additional control point in the computations, but requires no data. Section 2 will be assumed to 
be the same as section 3, unless the user enters a separate cross section for that location. 
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The WSPRO program first computes a natural profile (No bridge structure) by performing a 
standard step calculation from section 4 (exit section) to section 3F (full valley section); then 
from section 3F to section 1 (approach section). The program then computes the profile with the 
bridge in place. In this case the energy balance is performed from section 4 to 3 (Bridge 
opening); then from section 3 to 2 (upstream bridge opening); and finally from 2 to 1. The total 
energy equation from the exit to the approach section can be written as follows: 



2 ^ 





+ h 


1(1-4) 


(14) 


Where: h, == 

V, - 

h4 = 

V4 = 

K 


Water surface elevation at section 1. 
Velocity at section 1. 

Water surface elevation at section 4. 
Velocity at section 4. 

Energy losses from section 1 to 4. 


Energy losses from section 1 to 4 are equal to friction losses from 1 to 4 and an expansion loss 
from 3 to 4. The incremental losses are calculated as follows: 


From Section 4 to 3 

Losses from section 4 to section 3 are based on friction losses and an expansion loss. Friction 
losses are calculated using the geometric mean friction slope times the straight line distance 
between sections 3 and 4 (distance B from Figure 8 ). The following equation is used for friction 
losses from 3 to 4: 


TO-4) 


BQ^ 


(15) 


Where K 3 and K 4 are the total conveyance at sections 3 and 4 respectively. The expansion loss 
from section 3 to section 4 is computed by the following equation: 
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(16) 


K - ^ [2^4 - «4 - 


2gA 


Where a and P are energy and momentum correction factors for nonuniform flow, and P 4 are 
computed as follows: 


a 


4 


K^IA^ 


(17) 


K^IA 


ttj and p 3 are related to the bridge geometry and are defined as follows: 



(19) 



( 20 ) 


Where C is an empirical discharge coefficient for the bridge, which was developed as part of the 
Contracted Opening method by Kindswater, Carter, and Tracy (USGS, 1953). 


From Section 2 to 3 

Energy losses from section 2 to 3 (through the bridge) are based on friction losses only. The 
following equation is used for calculating the friction loss: 
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( 21 ) 


T(2-3) 


= L. 


(2-3) ) 


Where is the distance between the bridge sections, and K 3 is the conveyance at section 3 
inside the bridge. 


From Section 1 to 2 

Energy losses from section 1 to 2 are based on friction losses only. The equation for computing 
the friction loss is as follows: 


7(1-2) 




( 22 ) 


Where is the effective flow length in the approach reach, and Kj. is the minimum of the 
conveyances inside the bridge. The effective flow length is computed as the average length of 20 
equal conveyance stream tubes (FHWA, 1990). 

Pressure flow and Weir flow are modeled in a similar manner as in the HEC-RAS software. 


MODEL DEVELOPMENT 

Because this study is a comparison between models, it is very important that the geometric data 
and model parameters be developed in a consistent manner. The key pieces of information are: 
cross-section geometry and location; Manning’s n values; reach lengths (left overbank, main 
channel, and right overbank); contraction and expansion coefficients; bridge geometry; 
ineffective flow areas; and the location of the main channel bank stations. Each of the models 
were developed using the same data where possible (i.e., the same cross sections and locations, 
the same n values, bridge geometry, etc.). The only differences between the models were the 
required coefficients that are used in solving the hydraulic equations in the vicinity of the bridge. 
This will be discussed in further detail under the description of model calibration. 

CROSS SECTION DATA AND LOCATIONS 

At all of the bridge sites, the corresponding hydrologic atlas provided cross section data that was 
surveyed along the stream, both sufficiently upstream and downstream from the bridge sites. 
Each data set included no fewer than eight cross sections. Additional cross sections were also 
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surveyed just downstream of the bridge and inside the bridge. In general, the surveyed cross 
sections were used at the locations provided. In addition to the surveyed cross sections, it was 
necessary to locate cross sections in the vicinity of the bridge, as required by the various models. 
These cross sections consisted of: an exit section(section 1 from figure 1); a cross section 
upstream of the bridge (section 3 from figure 1); and an approach section (section 4 from figure 
1). The cross sections were formulated by making a copy of the closest surveyed section, and 
then making any minor adjustments to the station and elevation points if needed. 

There is some question about the sharpness of the contraction and re-expansion of effective flow 
in the immediate vicinity of a bridge, and indeed this concept is an integral part of the bridge 
modeling methodologies. For this study, almost all of the models are based on the assumption 
that both the contraction and expansion reaches have a length equal to the width of the bridge 
waterway opening. Because this assumption is a requirement in the WSPRO model, the 
approach and exit cross sections were placed at the same locations in HEC-RAS and HEC-2. 
This was done in order to have consistency between the models. In addition, for these specific 
sites the assumption of one bridge opening width for the expansion and contraction reach lengths 
appears to give reasonable results. However, for three of the events, it was found that improved 
results could be achieved by moving the exit section to a greater distance downstream from the 
bridge. This was done for both HEC-RAS and HEC-2. 


MANNING’S N VALUES 

The USGS hydrologic atlases contain color coded land use maps for the area of the floodplain. 
Each color is assigned a range of values for the Manning roughness coefficient. The maps also 
recommend using the high end of the range for depths below 0.6 meters, the low end of the range 
for flow depths above 1.0 meter, and linear interpolation between those two depths. 

The approach for selecting n values was to have only as many different n values in the model as 
there are roughness types on the maps. The only exception to this is where the channel and the 
overbanks have been coded with the same roughness type, and the channel flow depth is much 
greater than the overbank depth. The information in the atlases was represented as much as 
possible in the models, with respect to where changes in the n values can be expected to occur. 


CHANNEL BANK STATIONS AND REACH LENGTHS 

In some of the study reaches, the atlases delineate the low-flow channel on the planimetric view 
of the atlas. The models for these reaches have the section bank stations set to correspond with 
the atlas low-flow channel, even at sections where this alleged channel does not include the 
lowest point in the cross section. For those sites that did not indicate a low-flow channel location, 
bank stations have been set at each section generally to include the sharpest and deepest incision 
in the channel, except where this incision appears to be a tributary channel. In many of these 
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cases, the bank station settings are inconsequential since a uniform n value has been used across 
the section. 

Channel reach lengths for the models were scaled from the atlases. The main channel reach 
lengths were measured along the invert of the main channel. The overbank reach lengths were 
measured along the perceived path of the center of mass of the overbank flow. In general, the 
center of mass of the overbanks was estimated as being 1/3 the distance from the main channel 
bank station to the edge of the water surface boundary. 


CROSS SECTION EFFECTIVE FLOW AREA 

For the sections immediately adjacent to bridges and those which incorporate an apparent area of 
'dead water,' the effective flow area at the cross section has been constrained within the 
floodplain. In most cases this has been done in a manner which accomplished the constraint 
without adding wetted perimeter at the effective flow boundaries, using the assumption that a 
water-to-water boundary induces no significant friction loss. Two methods can be used which 
have this effect: using extremely high n values in the ineffective areas (n = 100); or using the 
program’s ineffective flow option. 

BRIDGE CROSSING GEOMETRY 

Detailed information about the bridge (road embankment, bridge deck, abutments, and piers) was 
well documented in each of the USGS hydrologic atlases. Each atlas contained a surveyed cross 
section inside the bridge, as well as detailed drawings of the bridge opening. This information 
was used consistently between all of the computer programs. 


MODEL CALIBRATION 

At all of the study sites, reach geometry, discharge, and observed water surface data are given far 
enough downstream that a profile starting point can be selected which is outside of the influence 
of the bridge. The main calibration parameter in all of the models was Manning’s n values. 
Contraction and expansion coefficients were also estimated. 

MANNING’S N VALUES 

Calibration of Manning’s n values was done using the entire observed extent of each reach. 
Exceptions to this are reaches where there was a significant inaccuracy in the downstream 
portions of the model, which could not be rectified by reasonable calibration. It was suspected 
that some important changes in roughness existed which were not documented on the atlas. In 
such cases, the starting point was moved upstream beyond the troublesome spot. 


24 



The calibration of n values was done by raising or lowering the n values of the entire reach. 
There are some locations in the models where a change in the n value at a specific cross section 
would produce a much better fit to the observed data, but no indication is given in the atlas that 
the roughness has changed for that cross section. In these cases, we did not change the n value. 
In other words, we have avoided a section-by-section calibration of n values, since to calibrate n 
values in this way would produce a model where all other parameter effects are masked. All 
changes in n values, within a particular model, were constrained to the n value ranges provided 
by the atlases. Also, the same n values were used in each of the three models for the same data 
set. 

All the models involved in the study incorporate the set of n values which, following the basic 
rules mentioned above, provide the best fit to the observed data over the entire reach, with 
greater weight given to those portions of the reach not likely to be affected by the bridge (cross 
sections downstream of the exit section). The final n values were set after we were confident 
that the reach geometry was modeled reasonably based on the atlases. 


CONTRACTION AND EXPANSION COEFFICIENTS 

In general, between cross sections contraction and expansion coefficients were set to 0.1 and 0.3 
respectively. There is some uncertainty about the proper values for these coefficients in the 
vicinity of a bridge. At most of the sites in our study, the velocity heads are low, even at the 
bridge. Therefore the setting of these parameters had little effect on the solution. We have 
considered two possible discreet settings: 

1. Contraction Coefficient = 0.3 and Expansion Coefficient = 0.5, or 

2. Contraction Coefficient = 0.5 and Expansion Coefficient = 0.8 

For most of the bridge sites, values of 0.3 and 0.5 were used for the contraction and expansion 
coefficients, respectively. In general, the setting of the contraction and expansion coefficients 
did not influence the results to any significant degree, for these specific data sets. 


FINAL MODEL ADJUSTMENTS 

Once the calibration of n values and contraction and expansion coefficients was completed, the 
performance evaluation and model comparison was carried out by selecting the downstream 
starting cross section as that section which is nearest to, yet outside of, the assumed expansion 
reach downstream of the bridge. The starting water surface was set to the observed value at this 
location in order to ensure that any differences between computed and observed water surface 
elevations in the vicinity of the bridge were due solely to the bridge computations. For most of 
the models the downstream starting section was the exit section. In a few of the models it was 
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necessary to set the starting water surface elevation at the cross section just downstream of the 
exit section. This was done, for these sites, because the observed water surface information in 
the vicinity of the exit section was inconclusive. Therefore, it was deemed that it would be better 
to start at the next downstream section with a better known water surface elevation. 


MODELING RESULTS 

During the course of this study it was found that for HEC-RAS and HEC-2, the energy based 
bridge solution methods (Normal Bridge for HEC-2) consistently produced better results than 
any of the other methods available in these models. This is consistent with the fact that friction 
losses are the predominant factor in the change in energy through the bridges at these specific 
locations. Because of this, the results shown for HEC-RAS and HEC-2 are only for the energy 
based solution methods. 

The results from this study are presented in the form of tables and graphics. Table 1 shows 
summary results for all of the events at all of the bridge sites modeled. This table shows 
observed and computed water surface elevations at the approach section (section 4 from Figure 
1 ) of each bridge. 

Table 2 shows the average absolute error in water surface elevations based on a comparison of 
three locations at each bridge site. These locations were: just downstream of the bridge; at the 
approach section; and the next cross section upstream of the approach section. The average 
absolute error at each location is computed by taking the sum of the absolute value of the 
computed minus observed water surface at the three locations, then dividing by three. 

In addition to Table 1 and Table 2, Appendix B includes a separate table and graphic for each of 
the bridge sites. The tables include observed and computed water surface elevations at several 
locations. An example of this output is shown in Figure 9 and Table 3. Figure 3 shows the 
observed and computed water surface profile results for the Poley Creek bridge site. Table 3 has 
the corresponding numerical results. 
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Table 1 

Computed and Observed Water Surface Elevations at the Approach Section 


STUDY 

LOCATION 

Flow 

(cfs) 

Observed 

WSEL 

(ft) 

HEC-RAS 

HEC-2 

WSPRO 

WSEL 

(ft) 

Error 

(ft) 

WSEL 

(ft) 

Error 

(ft) 

WSEL 

(ft) 

Error 

(ft) 

Alexander Cr. 

5508 

88.4 

88.2 

-0.2 

88.1 

-0.3 

88.3 

-0.1 

Alexander Cr. 

9500 

90.2 

90.1 

-0.1 

90.0 

-0.2 

90.1 

-0.1 

Beaver Cr. 

14000 

217.8 

217.9 

0.1 

217.8 

0.0 

217.3 

-0.5 

Bogue Chitto 

25000 

337.3 

337.8 

0.5 

337.5 

0.2 

337.6 

0.3 

Bogue Chitto 

31500 

338.3 

338.9 

0.6 

338.5 

0.2 

338.8 

0.5 

Buckhom Cr. 

4150 

322.0 

322.1 

0.1 

322.2 

0.0 

322.3 

0.3 

Cypress Cr. 

1500 

116.1 

115.8 

-0.3 

115.7 

-0.4 

115.9 

-0.2 

Flagon Bayou 

4730 

76.3 

76.2 

-0.1 

76.2 

-0.1 

76.9 

0.6 

Okatama Cr, near 
Magee 

16100 

367.2 

367.3 

0.1 

367.1 

-0.1 

367.3 

0.1 

Okatama Cr. East 
of Magee 

12100 

371.9 

371.5 

-0.4 

371.5 

-0.4 

372.6 

0.7 

Pea Cr. 

1780 

359.1 

358.9 

-0.2 

358.8 

-0.3 

359.4 

0.3 

Poly Cr. 

1900 

234.8 

234.7 

-0.1 

234.6 

-0.2 

235.0 

0.2 

Poly Cr. 

4600 

237.2 

237.2 

0.0 

237.2 

0.0 

237.6 

0.4 

Tenmile Cr. 

6400 

110.9 

111.0 

0.1 

111.0 

0.1 

110.9 

0.0 

Thompson Cr. 

3800 

200.3 

200.6 

0.3 

200.6 

0.3 

200.9 

0.6 

Yellow River 

2000 

234.2 

234.2 

0.0 

234.1 

-0.1 

234.3 

0.1 

Yellow River 

6603 

237.3 

237.7 

0.4 

237.5 

0.2 

237.8 

0.5 


Average Absolute Error 0.21 0.18 0.32 


27 


































































































































































Table 2 

Average Absolute Error in Water Surface Elevation Based on Three Locations 


STUDY 

LOCATION 

Flow 

(cfs) 

HEC-RAS 

HEC-2 

WSPRO 

Rank 

Ave Abs 
Error(ft) 

Rank 

Ave Abs 
Error (ft) 

Rank 

Ave Abs 
Error (ft) 

Alexander Cr. 

5508 

2 

0.20 

3 

0.30 

1 

0.17 

Alexander Cr. 

9500 

1 

0.10 

2 

0.13 

3 

0.23 

Beaver Cr. 

14000 

2 

0.07 

1 

0.03 

3 

0.50 

Bogue Chitto 

25000 

1 

0.23 

3 

0.27 

1 

0.23 

Bogue Chitto 

31500 

3 

0.40 

1 

0.27 

2 

0.33 

Buckhom Cr. 

4150 

1 

0.20 

1 

0.20 

3 

0.23 

Cypress Cr. 

1500 

2 

0.23 

3 

0.30 

1 

0.17 

Flagon Bayou 

4730 

1 

0.27 

2 

0.30 

2 

0.30 

Okatama Cr. near 

Magee 

16100 

1 

0.17 

2 

0.20 

3 

0.23 

Okatama Cr. East of 
Magee 

12100 

1 

0.40 

2 

0.47 

3 

0.80 

Pea Cr. 

1780 

1 

0.37 

1 

0.37 

3 

0.50 

Poly Cr. 

1900 

1 

0.30 

2 

0.37 

3 

0.40 

Poly Cr. 

4600 

1 

0.20 

2 

0.23 

3 

0.50 

Tenmile Cr. 

6400 

2 

0.23 

2 

0.23 

1 

0.10 

Thompson Cr. 

3800 

1 

0.20 

1 

0.20 

3 

0.40 

Yellow River 

2000 

1 

0.17 

1 

0.17 

3 

0.23 

Yellow River 

6603 

2 

0.40 

1 

0.30 

3 

0.43 


Mean Average Absolute Error 

0.24 

0.26 

0.33 

Total Number of No. 1 Rankings 

11 

7 

3 

Total Number of No. 2 Rankings 

5 

7 

2 

Total Number of No. 3 Rankings 

1 

3 

12 
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Figure 9. Computed and Observed Water Surface Profiles For Poley Creek. 


Table 3 

Computed and Observed Water Surface Elevations for Poley Creek. 


Cross 

Section 

Number 

Channel 

distance 

(ft) 

Obs. 

WSEL 

(ft) 

HEC-RAS 

HEC-2 

WSPRO 

WSEL 

(ft) 

Error 

(ft) 

WSEL 

(ft) 

Error 

(ft) 

WSEL 

(ft) 

Error 

(ft) 

4 EX 

1985 

234.9 

234.9 

0.0 

234.9 

0.0 

234.9 

0.0 

4.1 

2187 

235.3 

235.9 

0.6 

235.9 

0.6 

236.1 

0.8 

BR 

2189 

N/A 

236.0 


235.9 


N/A 


BR 

2217 

N/A 

236.6 


236.6 


N/A 


Qmi 

2219 

237.0 

236.7 

-0.3 

236.7 

-0.3 

N/A 


SAP 

2421 

237.2 

237.2 

0.0 

237.2 

0.0 

237.6 

0.4 

6 

3501 

238.0 

238.0 

0.0 

237.9 

-0.1 

238.3 

0.3 

7 

4341 

239.0 

239.0 

0.0 

238.9 

-0.1 

239.2 

0.2 
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SUMMARY AND CONCLUSIONS 


The purpose of this study was to evaluate the effectiveness of the HEC-RAS bridge modeling 
approaches at sites with extensive observed data, and to compare HEC-RAS with HEC-2 and 
WSPRO with respect to bridge modeling performance. Data were collected from the USGS in 
the form of Hydrologic Atlases (USGS, 1978). Seventeen flood events were analyzed at 13 
different bridge sites (four sites had two flood events). Hydraulic models of the bridge site, as 
well as several cross sections upstream and downstream, were developed using HEC-RAS, HEC- 
2, and WSPRO. Results were presented in tabular and graphical forms. Conclusions from this 
study are as follows: 

1. In general, all models were able to compute water surface profiles within the tolerance of 
the observed data. The observed data at each cross section is based on a series of measured 
points as well as high water marks. The variation of the water surface at any given cross section 
was on the order of 0.1 to 0.3 feet. The observed water surface was taken as the mean value of 
the active flow area. The mean absolute error in computed versus observed water surface 
elevations varied from 0.24 (HEC-RAS) to 0.33 (WSPRO). Given the small variance in the 
results, it is concluded that any of the models can be used to compute adequate water surface 
profiles at bridge locations, and that no one model performed significantly better than another. 

2. The HEC-RAS and HEC-2 energy based bridge solution methodologies are very similar 
to the WSPRO bridge procedure. All of the programs solve the water surface profile through the 
bridge by performing an energy balance. The following is a list of the differences between HEC- 
RAS, HEC-2, and WSPRO: 

By default HEC-RAS and HEC-2 use the average conveyance method for 
computing friction losses, while WSPRO uses the geometric mean method. All of 
the programs allow the user to select from four possible methods. 

HEC-RAS and HEC-2 compute expansion losses as a coefficient times the 
absolute change in velocity head from the section just downstream of the bridge to 
the exit section. WSPRO uses an expansion loss equation that was derived from 
an approximate solution of the momentum, energy, and continuity equations. The 
derivation is based on a rectangular section expanding into a wider rectangular 
section. 

WSPRO uses the stream tube concept to compute a flow weighted reach length in 
the contraction reach, while both HEC-RAS and HEC-2 compute a flow weighted 
reach length based on the user entered left overbank, main channel, and right 
overbank reach lengths. These reach lengths are used in computing the friction 
losses in the contraction reach. 

Both HEC-RAS and HEC-2 compute a contraction loss upstream of the bridge as 
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a contraction coefficient times the change in velocity head from the approach 
section to the section just upstream of the bridge. WSPRO does not compute a 
contraction loss in this reach, it is assumed to be zero. 

3. The following factors have been found to be far more important than which model is 
used: accurate cross section information; the placement of and location of cross sections in the 
vicinity of the bridge; Manning’s n values; adequately modeling the bridge geometry; defining 
ineffective flow areas in the vicinity of the bridge; and the contraction and expansion 
coefficients. 

4. Most of the large errors for all three models occurred adjacent to the downstream and 
upstream faces of the bridge. In most cases, all of the models did well in predicting the water 
surface at the approach section and further upstream. The errors in predicting water surface 
elevations adjacent to the bridge could be due to many factors. One factor may be an inadequate 
representation of the true active flow area in the vicinity of the bridge opening. Placing the 
ineffective flow area limits too far in, or too far away, from the edges of the bridge opening, can 
have a large impact on computing the water surface elevation. Another factor is that the flow 
may be highly multidimensional rather than one-dimensional, as being modeled by all three of 
these programs. Another possibility is that the downstream water surface may be controlled 
more by the bridge geometry than downstream conditions (i.e., model does not know about 
bridge geometry when computing downstream water surface). Also, the transition of the active 
flow area in the expansion reach is not a straight line between the downstream section and the 
exit section, this can have a large effect on the prediction of friction losses in the expansion 
reach, thus greatly effecting the computed water surface just downstream of the bridge 

5. The data for this study were all of a similar type, in that they all had wide and heavily 
vegetated floodplains. Manning’s n values for the overbanks varied from 0.07 to 0.24. 

Velocities in the main channel were low, and the variation of velocity in the vicinity of the bridge 
was small. Because of these factors, the flow expansion transition was found to occur more 
rapidly and over a shorter distance than the traditional 4:1 expansion rate. Most of the data sets 
were found to expand at a rate of 1 times the average obstruction length of the bridge roadway 
embankment. A few of the data sets required a longer expansion reach. The use of a 4:1 
expansion ratio would have resulted in over predicting the water surface elevations in the vicinity 
of the bridge. In general, placing cross sections too far apart in the expansion reach will cause an 
over prediction in the friction losses, thus causing the over prediction of the water surface 
elevations. The program uses an average friction slope times a weighted reach length to compute 
the friction losses. When the cross sections are too far apart, the averaging normally ends up in 
over predicting the friction losses. In the case of the bridge expansion reach, the friction slope 
(slope of the energy grade line) just downstream of the bridge is very steep. The flow will 
normally expand very rapidly at first, then the expansion rate will decrease. This causes the 
greatest changes in the friction slope to occur near the bridge. Because the transition of flow and 
change in the friction slope are not linear, none of the friction slope averaging techniques can 
adequately model the friction losses over a long reach. Thus the placement of the exit cross 


31 


section is very important in modeling the friction losses that occur in the expansion reach. Given 
that a limited number of data sets were used in this study, and that they all had several similar 
characteristics, generalizations about the appropriate expansion ratio, for all cases, could not be 
drawn. 

6 . At most of the bridge sites the contraction and expansion coefficients did not play an 
important role in the calibration and modeling of the water surface profiles for these bridge sites. 
Because the change in velocities in and around the bridges were small, the magnitude of 
contraction and expansion losses were generally small compared to the friction losses. At a few 
of the bridge sites, the contraction and expansion losses were significant in comparison to the 
friction losses. At all of these sites, traditional values of 0.3 and 0.5 were used for the 
contraction and expansion coefficients respectfully. These coefficients were not used as a 
calibration parameter for this study. However, in many studies these coefficients could play a 
much greater role in predicting the energy losses and water surface elevations in the vicinity of 
the bridge. 


32 



APPENDIX A REFERENCES 


Bradley, J.N., 1978. Hydraulics of Bridge Waterways, Hydraulic Design Series No. 1, Federal 
Highway Administration, U.S. Department of Transportation, Second Edition, revised March 1978, 
Washington D.C. 

Eichert, B.S. and J.C. Peters, 1970. Computer Determination of Flow Through Bridges, Hydrologic 
Engineering Center, Technical Paper 20, U.S. Army Corps of Engineers, Davis, CA. 

FHWA, 1990. User’s Manual for WSPRO - A computer model for water surface profile 
computations, Federal Highway Administration, Publication No. FHWA-IP-89-027, 177 p. 

Hydrologic Engineering Center, 1991. HEC-2, Water Surface Profiles, User's Manual, U.S. Army 
Corps of Engineers, Davis CA. 

Hydrologic Engineering Center, 1995. HEC-RAS, River Analysis System, User's Manual, U.S. Army 
Corps of Engineers, Davis CA. 

Hydrologic Engineering Center, 1995. HEC-RAS, River Analysis System, Hydraulic Reference 
Manual, U.S. Army Corps of Engineers, Davis CA. 

Koch-Carstanjen, 1962. Von de Brewegung des Wassers und Den Dabei Aufretenden Kraten, 
Hydrofynamim, Berlin. A partial translation appears in Appendix 1, “Report on Engineering Aspects 
of the Flood of March 1938," U.S. Army Engineering District, Los Angeles, May 1939. 

U.S. Army Corps of Engineers, 1963. Hydraulic Design of Reservoir Outlet Structures, EM 1110-2- 
1602, Washington D.C. 

uses, 1953. Computation of Peak Discharge at Contractions, Kindsvater, C.E., R.W. Carter, and 
H.J. Tracy, U.S. Geological Survey Circular 284, Washington, D.C., 35 pp. 

USGS, 1978. Hydrologic Investigation Atlases HA 591 - HA 611, Department of the Interior, 
Denver, CO. 

Yarnell, D.L., 1934. "Bridge Piers as Channel Obstructions," Technical Bulletin 442, U.S. 
Department of Agriculture, Washington D.C. 


A-1 



APPENDIX B 


Model Comparison Results for: 
HEC-RAS, HEC-2, and WSPRO 


This Appendix contains detailed tables and graphics for all of the bridge sites and flood 
events used in this study. Part 1 of this appendix contains the results from HEC-RAS (Energy 
Method), HEC-2 (Normal Bridge Method), and WSPRO. 

Part 2 of this appendix includes comparison tables and graphics for three events in which 
the observed water surface elevations were above the bridge low chords. This portion of the 
appendix was used to compare the alternative bridge modeling methods in HEC-RAS and HEC-2 
for high flows. This section contains results for the HEC-RAS energy method, HEC-RAS 
pressure and weir flow, HEC-2 Normal Bridge, and HEC-2 Special Bridge. The comparisons of 
high flows show no significant differences between the four methods for the three events that 
were analyzed. However, three data sets are far to few to make any judgments about the 
modeling methods for high flows. 
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Standard Error is equal to the square root of the Relative Error Squared divided by the number of RES values summed: (Rel. Error Sq / n)''0.5 
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Indicates that the observed water surface value is questionable either because there were very few values recorded on the atlas and/or because those that 
have been provided are not located within the active flow area. 

Relative Error Squared is equal to the square of the absolute error divided by the depth : (Abs. Error / Depth )^2. 

Standard Error is equal to the square root of the Relative Error Squared divided by the number of RES values summed: (Rel. Error Sq / n)^0.5 
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square root of the Relative Error Squared divided by the number of RES values summed: (Rel. Error Sq / n)'^0.5 
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Relative Error Squared is equal to the square or the absolute error divided by the depth : (Abs. Error / Depth )'^2. 

Standard Error is equal to the square root of the Relative Error Squared divided by the number of RES values summed; (Rel. Error Sq / n)''0.5 
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have been provided are not located within the active flow area. 

Relative Error Squared is equal to the square of the absolute error divided by the depth : (Abs. Error / Depth )^2. 

Standard Error is equal to the square root of the Relative Error Squared divided by the number of RES values summed: (Rel. Error Sq / n)'^0,5 
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Relative Error Squared is equal to the square of the absolute error divided by the depth : (Abs. Error / Depth )^2. 

Standard Error is equal to the square root of the Relative Error Squared divided by the number of RES values summed: (Rel. Error Sq / n)'^.5 
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Relative Error Squared is equal to the square of the absolute error divided by the depth : (Abs. Error / Depth ^2. 

Standard Eiror is equal to the square root of the Relative Error Squared divided by the number of RES values summed: (Rel. Error Scj / n)''0.5 
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Indicates that the observed water surface value is questionable either because there were very few values recorded on the atlas and/or because those that 
have been provided are not located within the active flow area. 

Relative Error Squared is equal to the square of the absolute error divided by the depth : (Abs. Error / Depth )'^2. 

Standard Error is equal to the square root of the Relative Error Squared divided by the number of RES values summed: (Rel. Error Sq / n)^0.5 
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Relative Error Squared is equal to the square of the absolute error divided by the depth : (Abs. Error / Depth )'^2. 

Standard Error is equal to the square root of the Relative Error Squared divided by the number of RES values summed: (Rel. Error Sq / n)''0.5 
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square root of the Relative Error Squared divided by the number of RES values summed: (Rel. Error Sq / n)^.5 
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Distance Along Channel^ ft. 
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Standard Error is equal to the square root of the Relative Error Squared divided by the number of RES values summed: (Rel. Error Sq / n)'^0.5 
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BRIDGE HYDRAULICS - USGS HYDROLOGIC ATLAS SITES 8/3/94 [TENM.XLW]CHART.XLC 



Distance Along Channel, ft. 
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BRIDGE HYDRAULICS > USGS HYDROLOGIC ATLAS SITES 6/13/95 [THOM.XLW]CHART.XLC 



Distance Along Channel, ft. 
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BRIDGE HYDRAULICS - USGS HYDROLOGIC ATLAS SITES 6/13/95 [YELLOW.XLW]CHART.XLC 



Distance Along Channel, ft. 





















BRIDGE HYDRAULICS - uses HYDROLOGIC ATLAS SITES 6/13/95 [YELL.XLW1ANALYSIS.XLS 






I - 
^ g 

O -p 






B 


3 

00 


CO 

wl 

o 

CQ 

s 

CO 

Pu 

O 

i 

! 

A- 

X 

PQ 


? 

'a 


’B 

tD 

'5 


I 


C -3 


§ 


*o 

o 

o 

c: 


•5 o 
2 •£ 


^ Id 


■5 

H, 


o. S 

X > 

« to 

g s 

S i3 
-C 1) 

2 B 
*o 2 

(O 

£ « 

3 §• 
a " 

to >» 

B 

I 

^ X 
^ 2 


<u ” 

•s ^ 

2 I 

I ^ 


3 3 


So -o 

•O <L> — 

•n 4) 


B 

o 

tfc 


« -^2 


§ c2 

S 2 
is 

? c -a 


O. « 

I! 


•o 

§ '-5 

♦S <U 
O 2 


2 2*- 


2 « 
£ 1 


C i 

O ^ 

w ^ 

a 5 


CO U 

S ;S 


w 

CO 

X) 

< 


2 ^2 
a, 5 


•o 

4> 


:2 ’cs 

i 


J> 2 
,£] 3 

s ^ 

o o 
'S tn 
^ >11 
(O 4> 

3 .2^ 
or w 


§ 

U 


S 5 


^ .S 

s S 
2 ^ 

« 2 


P 2 

I I 


c; '*-' 
3 o 

3 ^ 

CT O 
CO o 

di •« 

2 i 

-- ^ 


(U 

2 CO M 

(rt rs ^ 


O w 


c3 *0 « 


« X) 


Q Q 2 2 


2 i 

cS w 
> 1 
I 1 

0^ CO 


X 

w 


p!^ a. < 
CQ < Z 


B-34 






BRIDGE HYDRAULICS - USGS HYDROLOGIC ATLAS SITES 6/13/9S [YELL.XLW]CHART.XLC 



Distance Along Channel, ft. 
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CDMPARISaV: HEC Ifigh Flow Methods 

HEC-2 Normal Bridge 
HEC-2 Special Bridge 
RAS Energy Only 
RAS Pressure and Weir 
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BEAVER CREEK, LOUISIANA Q = 14000 cfs 



B-38 


have been provided are not located within the active flow area. 

Relative Error Squared is equal to the square of the absolute error divided by the depth : (Abs. Error / Depth )^2. 

Standard Error is equal to the square root of the Relative Error Squared divided by the number of RES values summed: (Rel. Error Sq / n)"'0.5 
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Distance Along Channel, ft. 
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Relative Error Squared is equal to the square of the absolute error divided by the depth : (Abs, Error / Depth )'^2. 

Standard Error is equal to the square root of the Relative Error Squared divided by the number of RES values summed; (Rel. Error Sq / n)^.5 
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Distance Along Channel, ft. 
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